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Abstract. We discuss the bulk properties of QGP produced at RHIC obtained 
at time of hadronization. We argue that hadronization of quark-gluon plasma 
occurs at a critical pressure near to 82 MeV/fm 3 , obtained for SPS energy range. 



1. Introduction 

We study the physical properties of the hot and dense fireball, and in particular, 
the thermal pressure at time of hadronization. In our approach, we consider that at 
RHIC, as well as at the top SPS energies, a local domain of thermally colliding gluons 
and quarks has been formed, the quark-gluon plasma (QGP). In the early stages of 
the collision a high density, and thus pressure buildup occurs which is followed by a 
fast relativistic collective matter outflow. This system expansion dilutes the density 
beyond phase equilibrium transformation to hadron gas (HG). 

This sub-dense system is unstable and can experience a sudden breakup, 
converting rapidly from quarks into hadrons, with free streaming particles, only 
'strong' hadron resonances are subject to possible further interactions. However, this 
does not alter the final stable particle yields. Thus, fitting the yields of particles 
using the statistical hadronization model (SHM), we can infer from the measured 
hadronic particle yields also the yields of other unmeasured hadrons. Summing the 
contributions of many gas fractions, we obtain the physical conditions of the fireball. 

In such a procedure, it is necessary to allow for greatest possible flexibility in 
characterization of particle phase space, consistent with conservation laws and related 
physical constraints of QGP breakup. In particular, the QGP yield of strange and 
light quark pairs has to be nearly preserved while QGP particles are distributed into 
final state hadrons. This is accomplished using parameters which describe the quark 
pair yields, i.e., j q and 7 S . While these parameters in both phases are not equal, i.e., 
1? GF 7^ 7i HG i ' — 1i s > the yield of pairs are similar, i.e., iV^ GP ~ N^ G , provided 
that the QGP breakup process is rapid, and that is seen in the HBT data. A value 
of 7i ^ 1 allows to control the density of particles at given temperature, but unlike 
the chemical potential, 7, acts in the same direction for particles and conjugate (anti) 
particles. How this works will become clear in next section. 

A jump-up in the phase space occupancy parameter 7^ replaces an increase in 
volume associated with a slow re-equilibrating hadronization with mixed phase, a 
reaction picture incompatible with many reaction observables, including HBT. The rise 
in occupancy, just like the rise in the volume size needed when chemical equilibrium 
7? G = 1 is assumed, accommodates transformation of a entropy dense QGP phase 
into entropy dilute HG phase. Similarly, a jump-up in the strangeness phase space 
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occupancy parameter 7 S allows for higher density of strangeness in QGP compared to 
HG. We have, in general, 

7 l HG (^/)>7 l QGP (i/), i = Q,s. 

The available number of quark pairs in QGP at hadronization decisively influences 
the possibility to form baryons. In fact, the ratio of baryon to meson yield arising 
in microscopic dynamics of hadronization is proportional to j q . This establishes 
the necessity to include the occupancy parameters in order to describe the yields 
of hadrons, since this is the parameter which allows for a hadronization dependent 
dynamical relative yield of mesons and baryons. Conversely, a study of particle yields 
with a fixed light quark equilibrium value 7^ G = 1 presumes that the relative yield 
of baryons to mesons is fully chemically equilibrated, and that we know well the 
spectrum of hadrons. Clearly, neither assumption is safe and the choice 7^ G = 1 is 
over-constraining any hadronization model, as of course is the choice 7^ G = 1. 

2. Particle yield SHM data analysis 

The analysis of experimental hadron yield results requires a significant book-keeping 
and fitting effort in order to allow for resonances, particle widths, full decay trees 
and isospin multiplet sub-states. A program SHARE (Statistical HAdronization 
with REsonances) suitable to perform this data analysis is available for public 
use [TJ [5] • This program implements the PDG [3] confirmed (4-star) set of particles 
and resonances, and we use [4] already for two years the modern cr-meson mass [5] 
{m a = 484, r CT /2 = 255 MeV). 

The important parameters of the SHM, which control the relative yields of 
particles, are the particle specific fugacity factors A and the space occupancy factors 
7 discussed above. The fugacity is related to chemical potential fx = Tin A. The 
occupancy 7 is, nearly, the ratio of produced particles to the number of particle 
expected in chemical equilibrium, and thus, meson yield is (nearly) proportional to 
7 2 and baryon yield to 7 3 (here, we did not distinguish the valance quark content for 
u, d, s quarks). The actual formula for the momentum distribution is, both for the 
HG and QGP phase: 

''" X m= J„ -n-,^,^ ^, (1) 



d 3 pd 3 x ~ J ^' (2tt) 3 7 - 1 A- 1 e £ / T ± 1 

where the Boltzmann limit of the Fermi '(+)' and Bose '(— )' distributions, applicable, 
in particular when m/T > 1, is indicated, g is the degeneracy factor, T is 
the temperature and E = E(p) is the single particle energy spectrum, typically 
E=^/m 2 +p 2 . 

The fugacity A is associated with a conserved quantum number, such as net- 
baryon number, net-strangeness or heavy flavor. Thus, antiparticles have inverse value 
of A, and A evolution during the reaction process is related to the changes in densities 
due to dynamics, such as expansion. Contrary to A, 7 is the same for particles and 
antiparticles. Its value changes as a function of time, even if the system does not 
expand, for it describes buildup in time of the particular particle species. For this 
reason, 7 changes rapidly during the reaction, while A is more constant. It is 7 which 
carries the information about the time history of the reaction and the precise condition 
of particle production referred to as chemical freeze-out. 
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Figure 1. x 2 /dof (top) and the associated significance level Pc[%] (bottom) as 
function of jq, on left for the AGS/SPS energy range, and on right for RHIC. 



In the quark phase, each particle has its proper chemical yield co-factor, thus for 
light quarks q — u, d, we have yield co-factors 7°" GP A°- GP , and for antiquarks q, we 
have 7 - P A~ 1C ^ , and similarly for strange quarks s and antiquarks s. In the HG 
phase, we need to count valance quark content of each hadron. For example, for the A, 
the chemical co-factor is jf G 7 2 HG A^ G A 2HG , while for A, it is 7? G 7g HG A;r 1HG A~ 2HG . 
We recall that the chemical potentials of baryon number, /zb, and hyperon number, 
MS, are 

/iB = 3T In X q , = T In \ q — T In A s , 

where, for historical reasons, hyperon number has opposite quantum number to 
strangeness. Above and from now on, when the upper index is absent the (chemical) 
variable considered refers to the final state phase, thus to HG. SHARE allows the 
conservation of (electrical) charge Q, which is done at the cost of introducing the 
fugacity A/3 = X u /^d and we note, in this context, that \ q = \J \ u \d- 

We evaluate the success of our data fit considering the profile of x 2 /dof as function 
of 7 9 . (see top of figure[T]), on left for AGS-SPS energy range, and on right for RHIC. 
The best fit is clearly not obtained at j q = 1 Note that for a small number of dof, the 
value of % 2 /dof can be very misleading, for this reason the bottom frames in figure [T] 
show the confidence level Pc = CL. 

The meaning of Pc = CL is explained in "Review of Particle Physics" [3] where 
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we see, in Figure 32.2, lines of fixed value of Pc, for given values x 2 /dof and dof. 
Note that the value Pc = 50% is equivalent to x 2 /dof = 1 for the case that very 
many dof are present. However, for a small number of dof, a very much smaller value 
of % 2 /dof must be achieved to claim good confidence fit. Only Pc(x 2 , dof), and not 
X 2 /dof, expresses confidence in the validity of the model used to fit the data. 

Effectively, Pc(x 2 , dof) also expresses confidence in the data, provided that we 
believe in the model. This is easily recognized by checking what happens when we 
intentionally alter an experimental data point, e.g., by 2 s.d.. We find that our data 
fit remains stable in the sense that we find nearly the same model parameters, but 
Pc becomes much smaller, and the falsified data point contributes dominantly to the 
error of the fit. 

At SPS (left side, bottom frame of figure [T]), our fits, carried out with the NA49 
2008 data set, converge to a common best value Pc = 70% [6j. This value has 
somewhat higher Pc than one should expect on statistical grounds (50%). However, 
we have treated systematic errors as if these were statistical errors, and have added 
these linearly to the statistical errors, in effect making many measurement errors too 
large. On the right hand side, in figure [TJ we see RHIC results. It seems that for 



^/snn = 200 and 130 GeV, the confidence level is here too high, suggesting that the 

combination of statistical errors with systematic errors was not appropriate (our, here 
presented, 62 GeV fit comprises some extrapolated and interpolated results and should 
not be seen as yet to be a 'real' data fit). 

3. Data used and statistical parameters 

The data sets for AGS-SPS study were presented elsewhere [BJ, the study of total 
(as opposed to dN/dy) hadron yields at RHIC is relying on extrapolations, and the 
discussion of this matter goes beyond the scope of this presentation. In table [IJ we 
show the in-out data fields for the fits to central rapidity 200, 130 and 62.4 GeV at 
RHIC, as used here. As seen in table [TJ we combine PHENIX data for direct single 
particle spectra with RHIC date for (strange) particles yields reconstructed from the 
decay particles using, e.g., the invariant mass method [7J, a more complete discussion 
of our data set goes beyond the scope of this report, see also [4]. 

Other statistical hadronization (SH) parameters we derive from the data, shown 
in bottom of table [TJ are the source volume V (that is dV/dy for RHIC, the 
volume associated with the interval of rapidity in which particles are measured), 
the temperature T, at which particles stop changing in yield (chemical freeze-out). 
Moreover, we obtain chemical potentials fiB = 3/i g = 3TlnA g , /ig = Tln(A g /A 5 ), 
related to conserved quantum numbers: baryon number and strangeness, respectively. 
We also obtain A/3 which expresses the asymmetry in the 3-rd component of the 
isospin. Especially for low energy reactions, where the particle yield is relatively low, 
this parameter differs significantly from unity. We have become aware by checking 
the work of other groups pursuing statistical hadronization of QGP and fits to hadron 
yields that the net charge per net baryon ratio (0.39 for heavy nuclei) is not maintained 
in this work. 

4. Hadronization condition 

The SHARE program provides, beyond statistical parameters, also an opportunity to 
evaluate the physical properties of the bulk matter at hadronization. These show a 
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Table 1. The constrains, imposed and natural at the top, followed by input 
particle data and the resulting statistical parameters, and at bottom, the chemical 
potentials derived from these, for RHIC central rapidity, most central collisions. 
A s values are obtained from the constraint to zero strangeness. The weak feeds 
allowed were as stated by experimental groups or/and estimated by us, however, 
in the fit of the p,p yields we accepted complete weak feed from all hyperons. 



V s nn [ijevj 


OZ.4 


1 QA 


OAA 

zUU 


P rfinV/fm 3 ! 

i [Vjre v / un j 


n aqo_i_a nni 


n aqo_i_a nm 


a aooia nm 


F,lh ffipVl 




a A7E 


a 7^4-n 07^ 


( q q\ /fa —1— a\ 


n 


n 


n 

u 


O ih 




0.4±0.01 


0.4±0.01 


7T+ 


233+26 


276+36 


286 4+24 2 


7V~ 


237±27 


270±36 


281.8+22.8 


K+ 


38±4.3 


46.7±8 


48.9+6.3 


K~ 


32.6±4.7 


40.5±7 


45.7+5.2 


d>/K" 




0.15±0.03 


0.174+0.03 


f 


34.3±3.8 


28.7±4 


28.3+4.8 


f 


13.8±1.6 


20.1±2.8 


13.5+1.8 


A 




17.35±0.8 


16.7+1.3 


A 




12.5±0.8 


12.7+1.1 




1.84±0.2 




2.17+0.25 




1 i 6+0 1 2 




1 83+0 25 


E~/h~ 




0.0077±0.0016 








0.853±0.1 






0.229±0.035 






0, 


0.176±0.030 






n + n 






0.85+0.08 


(Q + Q)/h~ 




0.0021±0.0008 




K (892)/K" 




0.26±0.08 


0.23+0.05 


dV/dy [fm 3 ] 


1089±74 


1172±93 


1156+88 


T [MeV] 


142.9±0.3 


140.2±0.2 


140.5+0.5 


\ 


1.166±0.036 


1.077±0.020 


1.066+0.030 


\s 


1.066* 


1.029* 


1.033* 


7<J 


1.53±0.10 


1.56±0.031 


1.54+0.14 


7s 


1.74±0.22 


2.32±0.31 


2.30+0.40 


A/3 


0.992±0.003 


0.997±0.001 


0.997+0.002 


fiB [MeV] 


65.9 


31.2 


27.1 


Ms [MeV] 


13.5 


6.4 


5.6 



change, from a low density and low pressure system at low y/s (AGS, lowest SPS 20 A 
GeV data) to a highly compressed phase just above this in energy. In figure[2j we show, 
in the upper frame, the pressure P we obtain for the different fits. We find [HE], in the 
study of the high energy SPS data, that hadronization is characterized by a remarkably 
constant value of P ~ 82 ± 2 MeV/fm 3 . This result arises in the SPS energy domain 
(which is involving the total particle yields) without a further constrain. For RHIC, 
the central rapidity results lack a fixed baryon number and, as shown in previous 
section, the errors are way too large yielding fit confidence which is too high. We thus 
decided, at RHIC, to introduce as a additional 'measurement' the value P ~ 82 ± 2 
MeV/fm 3 . 
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Figure 2. Pressure and energy per baryon as function of reaction energy. P is 
fixed for RHIC to the value at SPS to help find a reliable reproducible fit. 



In figure^ we also show in the lower frame the fireball energy per primary hadron, 
E/hp. This value is also remarkably constant for top SPS and all RHIC reaction 
energies. The variation we see is, in part, explained by baryon density variation, and 
at low energy by different properties of the hadronizing system also seen in many 
other observables. Note that if 7^ G = 1 is forced, the value E/hp — 0.75 rises 
to 1 GeV [8J. Considering the result achieved at SPS we introduce at RHIC as an 
additional constraint E/hp = 0.75 ± 10% as is also shown in table [T] With the two 
constraints, P and E/h p we find very good data fits with the outcome confirming that 
the hadronization pressure offers a good characterization of QGP breakup. 

Hadron particle pressure emerges, in our study, as a common physical property 
which defines when and how QGP breaks up into hadrons. Why this is so is explained 
remembering that, at P = 1 atm, water boils in New York and in Beijing at 100° C. 
The color non-conductivity of the true vacuum acts like a 'pot cover' keeping quarks 
together, the cover recedes when the pressure is high, QGP expands. After QGP 
breaks, the residual quark pressure turns into hadron pressure. In this picture 
the quark particle pressure has just the magnitude required to balance the vacuum 
pressure. Thus, the critical pressure of hadronization must be the vacuum pressure 
confining color. 

The pressure P is compared to several other physical bulk properties of QGP 
in figure [3l At low energy considering for example, entropy a, we note that the 
AGS and lowest SPS results agree and produce a low value, suggesting a source 
which is half as dense compared to other results. This would be just what one 
expects if QGP is not formed at low reaction energies, or/and when there has been a 
considerable re-equilibration of hadronization products. Consideration of other, more 
penetrating observables, such as strangeness per entropy s/S and the continuity of 
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Figure 3. Physical bulk properties at hadronization beginning with highest 
energy AGS up to RHIC 200 GeV results. Squares (blue) are total yield fits, 
triangle (red) are central rapidity data fits, all lines guide the eye, fit results are 
the symbols. 



total strangeness production as function of energy support the second hypothesis, a 
well equilibrating QGP fireball after hadronization. 

At high energy, we see that SPS and RHIC bulk properties are consistent. 
Moreover, the behavior as function of energy of, for example, the net baryon density 
(bottom frame) is consistent with the expectation that it should be decreasing — since 
baryon transparency increasing with energy of reaction is an intuitive requirement. 
Fits which force 7* = 1 can fail to produce this natural result. 

The values of energy density, E/V — > 500MeVfm~ 3 = (250 MeV) 4 , and entropy 
density, S/V — > 3.4 fm~ 3 , obtained for the high energy reactions are worth noting. 
These complement P -> 82 ± 2MeVfm~ 3 = (158 MeV) 4 . 

5. Particle yield predictions 

In figure[31 we show growth of strangeness pair yield with energy, squares (blue) for 
the total hadron most central reaction trigger, and RHIC central rapidity most central 
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trigger as triangles (red). In table [H we show the particle yields we obtain, comparison 
wit table Q] shows that our fit is very successful. 



Table 2. Output hadron multiplicity data for the RHIC energy range. See text 
for the meaning of predictions of yields at 62.4 and 130 GeV and of dN/dy at 
62.4 GeV. The input statistical parameters are seen in table [T] b = B — B = N~w 
for 47r results and b = d(B — B)/dN for results at central rapidity. Additional 
significant digits are presented for purposes of tests and verification. All yields 
are without the weak decay contributions. 
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6. Comments and Conclusions 

We believe that the assumption of j q = 1, we often see in literature in the context 
of the analysis of hadron particle yield data, tests the hypothesis that QGP was 
not, or only extremely briefly present in relativistic heavy ion collisions. The tacit 
assumption made is that instead of a QGP, a long-lasting cascade of hadronic reactions 
allows the HG particles to chemically equilibrate. However, if this were to be true, 
an even stronger evidence for HG dominance of HI reactions would be that the value 
7 g — > 1 emerges in the analysis, rather than being assumed. However, we find allowing 
7i 7^ 1, i = q, s that, nearly always, the chemical non-equilibrium prevails, with a 
much higher confidence level. Furthermore, this additional freedom produces QGP 
like properties of the bulk from which particles emerge. 

To summarize, we differ from other groups in the following aspects in our data 
analysis: 

(i) We use j q ^ 1, and thus allow the ratio of baryon to meson yields to be fixed 
independently of the hadronization temperature; 

(ii) We enforce, in the fit, the conserved ratio of charge to baryon number Q/b = 
0.395 ± 0.01, and are able to fit the associated A/3 fugacity; 

(iii) We do not enforce exact strangeness conservation, (s) — (s) = 0, but instead, we 
allow 5s = (s — s)/(s + s) = ± 0.05 to behave like a measurement, the reasoning 
is as follows: 

a) Summing all measured and unmeasured hadrons in strangeness 'conservation' 
condition, As = J2i K ~ J2j H ~^ 0) combines independent measurement errors 
and thus, even if the experimental data had all strangeness carrying hadrons, 
there would be a residual statistical error present in 5s; 

b) Some strangeness could escape detection in unknown 'particles', for example 
being bound in (nearly) uds-quark-symmetric semi-stable strangelett (a small 
drop of quark matter), this leads to 5s < — which is what we find as a preferred 
result in low energy fits; 

c) The experiments did measure many, but not all relevant particles carrying 
strangeness, e.g., S ± has not been measured, this yield is uncertain and thus 
5s = cannot be ever imposed on experimental grounds alone; 

d) Introducing an error in 5s, we perform a test of the hypothesis that weak decays 
remains weak in QGP phase, and thus, the net strangeness remains conserved — 
another way to understand this is to note that we cannot confirm that weak 
decays in QGP remain weak to better than the progressing error of individual 
contributing measurements. 

(iv) In the study of RHIC data, we consider a mix of STAR and PHENIX particle 
yield results, following the principle that the PHENIX integrated single particle 
spectra are more reliable than those of STAR, and STAR, in turn, given its high 
acceptance is more reliable in evaluation of yields of hadrons observed in their 
two or more particle decay channels. 

The reader should note that we have, in our fits compared to many other efforts, 
three more parameters: j q as explicitly stated above, A s since we do not fix but fit 
strangeness conservation and A/ 3 since, as the only group, we enforce a fit to Q/b. 
Thus, we have, in general, one less degrees of freedom (three more parameters and 
two more 'data' points, 5s = ± 0.05 and Q/b = 0.395 ± 0.01). However, the real 
issue is that all told we have a 7-dimcnsional space of parameters, X 1 , V,A^, A s , 'jq 7 'y s 
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and A/3, which contains many false minima, and the art of finding the domain of the 
best fit minimum is not easily acquired, and cannot be dispensed with the comment 
'fit is unstable'. Naturally, it takes much more effort to find a true minimum in a 
7-d parameter world, compared to a 2-d parameter world which corresponds to the 
simplest and least physical "equilibrium" model with T, /xb as parameters, i.e., setting 
also 7 S = 1, fixing /is by strangeness conservation and using only particle ratios in the 
fit. Whoever practices this today has learned nothing from the work of past 20 years. 

The most intriguing result of this analysis is the smoothness, and even near 
constancy, of physical properties of the fireball at chemical freeze-out condition seen 
for the top three SPS energies 40, 80 and 158 AGeV, which result agrees well with 
RHIC fits where we impose the pressure P and hadronization particle energy E/h p . Of 
particular physical interest is the value of hadronization pressure, P ~ 82 MeV/fm 3 , 
obtained at SPS, and found in this work to be consistent with RHIC data — one 
may imagine that in a phase transformation from quarks to hadrons, the pressure 
of quarks is transferred into the pressure of color-neutral hadrons, which can escape 
from the deconfined fireball. Since the flow pressure of quarks transfers smoothly into 
that of hadrons, we conclude that the thermal pressure of produced hadrons, P ~ 82 
MeV/fm 3 , provides a first estimate of the pressure of the vacuum which keeps color 
charged quarks inside the fireball up to the point of sudden fireball break-up. 

In summary, we presented a high confidence fit of high ccntrality data form AGS, 
SPS and RHIC, and have found common ground of results at SPS and RHIC regarding 
the bulk properties of hadronizing matter. This suggests that a deconfined phase is 
with great probability already formed at or near 30 A GeV. 

Acknowledgments 

Supported by a grant from the U.S. Department of Energy, DE-FG02-04ER41318 . 
Laboratoire de Physique Theorique et Hautes Energies, LPTHE, at University Paris 
6 and 7 is supported by CNRS as Unite Mixte de Recherche, UMR7589. 

References 

[1] G. Torrieri, S. Steinke, W. Broniowski, W. Florkowski, J. Letessier and J. Rafolski, Comput. 

Phys. Commun. 167 (2005) 229 [arXiv:nucl-th/0404083] . 
[2] G. Torrieri, S. Jeon, J. Letessier and J. Rafelski, Comput. Phys. Commun. 175 (2006) 635 

[arXiv: nucl-th /0603026] . 

[3] W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (20 06) 1. 

[4] J. Letessier and J. Rafelski, Eur. Phys. J. A 35 (2008) 221 arXiv:nucl-th/0504028] . 
[5] R. Garcia-Martin, J. R Pelaez and F. J. Yndurain, Phys. Rev. D 76 (2007) 074034 
arXiv:hep-ph/0701025 ; R. Kaminski, R. Garcia-Martin, P. Grynkiewicz and J. R. Pelaez, 

"Sigma pole position and errors of a once and twice subtracted dispersive analysis of pi-pi 

scattering data," arXiv:0811.4510 [hep-ph].In the proceedings of 14th High-Energy Physics 

International Conference in Quantum Chromodynamics (QCD 08), Montpellier, France, 7-12 

Jul 2008. 

[6] J. Rafelski and J. Letessier, arXiv:0901.2406 [hep-ph], in proceedings of Confinement 2008, 

Mainz, Germany, PoS(Confinement8)lll, (2009) M. Neubert, editor. 
[7] Anthony R Timmins, for the STAR Collaboration "Overview of Strangeness Production at the 

STAR Experiment" larXiv:0812.4080H [nucl-ex], in this volume. 
[8] J. Cleymans, H. Oeschle r, K. Redlich and S. Wheaton, Phys. Rev. C 73 (2006) 034905 

[arXiv:hep-ph/0511094] . 



